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We investigate the changes in the electronic structure of copper phthalocyanine (CuPc) crystals 
that is caused by intercalation with potassium. This is done by means of ab initio LSDA and 
LSDA+U calculations of the electronic structure of these molecular crystals. Pristine CuPc is found 
to be an insulator with local magnetic moments and a Pc-derived valence band with a width of 0.32 
eV. In the intercalated compound K2CuPc the additional electrons that are introduced by potassium 
are fully transferred to the Cg states of the Pc-ring. A molecular low spin state results, preserving, 
however, the local magnetic moment on the copper ions. The degeneracy of the Cg levels is split by 
a crystal field that quenches the orbital degeneracy and gives rise to a band splitting of 110 meV. 
Molecular electronic Coulomb interactions enhance this splitting in K2CuPc to a charge gap of 1.4 
eV. The bandwidth of the conduction band is 0.56 eV, which is surprisingly large for a molecular 
solid. This is line with the experimentally observation that the system with additional potassium 
doping, K2. 75 CuPc, is a metal as the unusually large bandwidth combined with the substantial 
carrier concentration acts against localization and polaron formation, while strongly promoting the 
delocalization of the charge carriers. 



Introduction Molecular crystals such as for instance 
fuUerenes, polyacetylene or (BEDT-TTF)2X salts which 
are chemically doped show a dramatic change in their 
electron transport properties Q • Starting out as insu- 
lators, upon doping they can become conducting and in 
special cases even superconducting. Because of the chem- 
ical richness of molecular crystals, building new metals 
by doping them remains an exciting and ever moving 
front. Very recently the Delft group showed that tran- 
sition metal phthalocyanines (MPcs) FePc, CoPc, NiPc, 
CuPc, initially insulating, can be turned into a metal 
through potassium doping 3]. These molecules are par- 
ticularly interesting because they have a magnetic mo- 
ment, which opens the way to interesting electron corre- 
lation effects. Because the molecules can have in addi- 
tion a possible orbital degeneracy, Jahn- Teller coupling, 
Hund's rule exchange and large on-site Coulomb inter- 
actions, the doped metallic molecular solids constitute a 
novel class of strongly correlated metals 0| . The question 
arises how upon doping the interplay develops between 
metallicity, Jahn- Teller distortions and collective correla- 
tion effects such as magnetic ordering. We set the stage 
for more involved many-body calculations by perform- 
ing ab initio electronic structure calculations on pristine 
and potassium doped copper phthalocyanine K2CUPC, 
which has the structure of the experimentally observed 
metallic phase with x — 2.75. For this doped compound 
K2.75CUPC we wish to establish: (i) whether the potas- 
sium atoms in donate their electrons to the CuPc, and if 
so, how many and in which molecular orbitals they oc- 
cupy (Cu-related or Pc-related), (ii) which are the prop- 
erties the valence and conduction bands that are de- 
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TABLE I: Electronic ground state configurations for CuPc 
with various valencies. The ai„ and Cg levels are Pc-derived 
and the big level is a Cu orbital. 



rived from these molecular orbitals, and in particular, 
(iii) whether the bandwidth such that it explains the ob- 
served metallic behavior of our doped molecular crystal. 
This indeed turns out to be the case: the metallicity of 
K2.75CUPC is due to its relative large bandwidth. 

Pristine CuPC: LSDA The structure of the CuPc 
molecule and its electronic energy level diagram |E 0| is 
shown in Fig. ^ The molecule is flat and has a four-fold 
rotation symmetry, D^h- CuPc is an open shell molecule: 
the Pc-ring is essentially a closed shell system, but copper 
is divalent and thus in a cP configuration, with a single 
hole that carries an uncompensated magnetic moment. It 
is important to note that the first electron removal state 
of the CuPc molecule is a Pc-derived aiu state and the 
first electron addition state is a Pc-derived Cg state. Also 
adding a second (and even third) added electron will not 
fill the Cu big state, but instead be added to the Cg state 
of the Pc-ring, see Table □ 

We start by investigating the electronic structure of un- 
doped solid CuPc before studying the effects of doping. 



FIG. 1: The CuPc molecule with symmetry. It is planar 
and the central Cu^^ ion is surrounded by two rings of each 
4 nitrogen atoms. Carbon and hydrogen make up the other 
atoms in the molecule. 



In its /3-structure the unit cell has a C2/1 space-group sym- 
metry and there are two CuPc molecules per the unit cell, 
the first one in the origin of the unit cell while the second 
molecule at (a/2,6/2,0) 8]. The molecules are slightly 
distorted and their rotation symmetry in the crystal is 
lowered to a two- fold rotation D2h- 

The electronic structure calculations are done within 
the density functional approach in the local density ap- 
proximation (LDA) and generalized gradient approxima- 
tion (GGA) for the intercalated system, using the Vienna 
Ab-initio Simulation Package (VASP) |3|. We first deter- 
mined the magnetic groundstate of the compound. From 
the computations in the SGGA framework we find that 
the ground state exhibits antiferromagnetic order, and 
that the energy gain of the antiferromagnetic state with 
respect to the paramagnetic one is 394 meV per unit cell. 
The magnetic moment of each molecule is close to 1 fis- 
However, we also find that the energy difference between 
the ferro and antiferromagnetic configuration of the two 
Cu spins in the unit cell is very small, i.e. smaller than 1 
meV. This indicates that even though local moments are 
present in this system, a magnetic ordering temperature 
-if ordering occurs at all- is expected to be very low. It 
also implies that the band structure that we will deter- 
mine next is nearly identical for the ferro and antiferro- 
magnetically ordered groundstate. In these total energy 
calculations a plane-basis set with a kinetic energy cutoff 
of 550 eV was chosen and the tetrahedron method with 
Blochl correction using 40 irreducibile k-points was used. 

The band structure and density of states of CuPc, as 
calculated by LSDA is shown in Fig. O The bands are 
arranged in almost degenerate pairs, indicating that the 
local electronic structure of the two molecules in the unit 
cell is practically identical. Moreover the bands show 
little dispersion, from which we conclude that the inter- 
action between the molecules is weak. In fact there is 
only appreciable dispersion of the bands along the crys- 
talographic 6-direction, as can clearly be seen from, for 
instance, the T — Y cut through the BriouUin-zone in 
Fig. 13 Thus from an electronic point of view, solid CuPc 
is quasi one-dimensional. 

In Fig. O we see that below the Fermi level there is 



FIG. 2: Schematic representation of the crystal structure of 
pristine CuPc. There are two inequivalent molecules in the 
unit cell. 



a very flat band, which is derived from localized Cu big 
states. Above the Fermi-level is the corresponding empty 
Cu state. At still higher energy are the empty Cg states. 
It is interesting to see that in the band structure the 
Cg states, which are doubly degenerate in the molecule, 
split up into two non-degenerate bands. The degeneracy 
is lifted in the crystal because of the ionic and covalent 
molecular crystal fields that distort the CuPc molecule. 

In contrast to the band calculation stands the molec- 
ular energy level scheme, from which we know that the 
first electron addition and removal states are not the Cu 
big states, but the Pc derived Cg and ai„ states, respec- 
tively. It is very unlikely in the solid this level scheme 
is altered very much. Rather, the different results in the 
solid are due to the fact that LSDA treats electron corre- 
lations between localized electrons insufficiently. Here we 
are dealing with molecules that have an open d shell and 
electron-electron correlation effects are known to be very 
strong for the Cu dP electrons. To include these effects 
one needs to go beyond LSDA. 

Pristine CuPc: LSDA + UWe incorporate the electron 
correlations that are due to the copper ions by an effective 
Coulomb interaction Ucu in the LSDA + U calculation 
scheme 0. The parameter Ucu is determined indepen- 
dently by computing the variation with respect to occu- 
pation number of the energy eigenvalue of the big state 
for a single relaxed CuPc molecule, using the GAMESS 
code [13. We find using a TZV basis set and the B3LYP 
functional that Ucu = 6.28 eV, a value that is typical for 
Cu . The presence of such strong interaction effects is 
in accordance with the recent calculations of Tosatti et 
al. '4^, which emphasize the importance of strong elec- 
tron correlations. Note that the parameter Ucu that we 
calculated includes screening of the Coulomb interaction 
by the electrons on the Pc-ring of the molecule. In a 
crystal there will be additional screening of Ucu by the 
surrounding molecules. This will effectively reduce Ucu 
[ill IT^ , but we find that our conclusions on the sym- 
metry, character and width of the valence and conduction 
do not depend on such additional screening. 

The LSDA+U bandstructure is shown in Fig. ??. The 



3 




FIG. 3: Band-structure, density of states and projected den- 
sity of states for CuPc, calculated with SGGA. The energies 
are plotted along lines in the first Brillouin zone, connecting 
the points: The Fermi-level is the zero of energy of the sys- 
tem. A = (0.5, 0.5, 0.0), F = (0.0, 0.0, 0.0), X = (0.5, 0.0, 0.0), 
L = (0.5,0.5,0.5), Y = (0.0,0.5,0.0), M = (0.0,0.5,0.5), 
Z = (0.0,0.0,0.5). 



system is magnetic, with an energy gain with respect to 
the paramagnetic state of 0.622 eV per molecule. Due 
to Ucu the flat 5ig-derived Cu band below the Fermi- 
level is shifted downward in energy with respect to the 
Pc derived states and the empty big is shifted upward. 
Now both the HOMO and LUMO derived band is of Pc 
character, in accordance with the energy level scheme of 
the bare molecule. We find a valence band width of 0.32 
eV. The correct order of the bands is very important, 
because next we will consider the electron doped system, 
for which the character of the valence band (very flat Cu- 
derived or Pc-derived and dispersing) is essential. From 
the calculations we find a bandgap of A^and = 1-2 eV. 
This is should be contrasted with molecular CuPc, which 
has an energy gap of AcuPc = IcuPc - AcuPc = 4.5 ± 0.3 
eV 0, 01, where / and A are the ionization potential 
and electron affinity, respectively. These two numbers, 
however, should only be compared when two additional 
effects are taken into account. First, due to the polariza- 
tion of the surrounding molecules the gap in the solid is 
screened with respect to its molecular value 0, 0, 0| . 
We calculated this screening energy to be 1.7±0.1 eV, so 
that on the basis of this result one would expect to find 
a bandgap in the solid of about A « 2.8 ± 0.3 eV. Thus, 
LSDA+U underestimates the gap. It is obvious that this 
is due to the fact that in such a calculation only the elec- 
tronic correlations for the Cu-electrons are taken into ac- 
count. The electron-electron interactions of the electrons 
on the Pc-ring, which makes out the largest part of the 
molecule, are not properly treated. We therefore deter- 
mined with GAMESS the Coulomb interaction, the Hub- 
bard t/pc, for electrons on the Pc-ring of the molecule. 
We find a value of Upc = 3.0 eV. Taking this into account 
we find a charge gap on the basis of the band calculations 



FIG. 4: Band-structure, density of states and projected den- 
sity of states for CuPc, calculated with SGGA-I-U. 



of A « 1.2 + 3.0 - 1.7 = 2.5 ± 0.1 eV. Thus within the 
error bars the two estimates of the gap agree and they 
are also consistent with the reported experimental val- 
ues from photoemission spectroscopy (2 3 ± 0.4 eV [T^). 
transport properties (2.1 ±0.1 eV |T3) and the lower 
bound on the gap of 2.1 eV ^B^l that is determined by 
energy electron loss spectroscopy. 

Doped Cu-phthalocyanine We used the structural re- 
finement of potassium intercalated CuPc reveals that the 
K2.75CUPC forms a stable phase. Its unit cell is mono- 
clinic, has C2h space-group symmetry, and contains two 
CuPc molecules ^ . The structure is quite different from 
pristine CuPc, because now there is one molecule in the 
unit cell at (0,0,0), as before, but the other one in the 
center of the cell, at (a/2, 6/2, c/2). According to the 
X-ray data there are two inequivalent positions of the 
potassium ions, which are both partially occupied. To 
perform calculations, we initially occupy only one of the 
two equivalent K-positions -the so-called K29- so that we 
have four potassium ions per unit cell, i.e. we start by an- 
alyzing K2CUPC before looking at the electronic proper- 
ties of K2.75CUPC. There are two reasons to occupy only 
one K-position. First, when investigating the electronic 
structure of our material we want to separate its intrinsic 
electronic properties from the ones induced by disorder 
on the K-sites. Second, the present unit cell is already 
very large (118 atoms) and the additional inclusion of 
disorder (for instance in a super-cell or coherent poten- 
tial approximation method) is numerically too demand- 
ing. We will see the band-structure of K2CUPC forms the 
proper starting point to understand the electronic struc- 
ture of the higher doped K2.75CUPC compound. 

Before being able to calculate the full bandstructure 
of K2CUPC , we need to determine its groundstate prop- 
erties. To this end, we first optimize the total energy 
allowing the ionic positions to relax. We observe that 
the K-ions move slightly in the unit cell to a position 
that is directly above the center of a phenyl ring of the 
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FIG. 5: Schematic representation of the crystal structure of 
K2.75CUPC. 



phthalocyanine molecules -a high symmetry position. In 
addition, we have also looked at the magnetic structure 
of the ground state. Again, these calculations are done 
with a kinetic cut-off energy of 550 eV and the tetrahe- 
dron with Bloch correction, using 40 irreducible k-points. 
The results of the SGGA + U band calculation on anti- 
ferromagnetic K2CUPC and the DOS are shown in Fig.El 

We see that the DOS of the valence and conduction 
bands does not have any K-character. We thus conclude 
that each potassium atom fully donates an electron to the 
CuPc rings so that we are dealing with CuPc^^. Apart 
from this, the K's do not play any role in determining the 
electronic properties. This implies that for a qualitative 
understanding of the electronic structure of potassium 
doped CuPc, the electronic structure of a CuPc molecule 
is a reasonable starting point. 

We focus in Fig. El^n the five bands around the Fermi 
level: the two that make up the valence band below it 
and the three of the conduction band above. Compared 
to pristine CuPc the valence band has taken up two ad- 
ditional electrons per molecule. This valence band is en- 
tirely derived from molecular eg orbitals, states that are 
localized on the Pc-ring. The eg orbitals of a neutral 
CuPc molecule can absorb up to four electrons because 
these states are two-fold degenerate (besides spin degen- 
erate). So for K2CUPC each CuPc molecule is left with 
two holes in the eg orbital. These empty states form two 
out of the three conduction bands. The third conduction 
band is the derived from the empty molecular big states, 
which is centered on the copper ions and has very little 
dispersion. 

From Fig. we see that the valence and conduction 
band, which are derived from the molecular eg states, are 
separated by a small gap of 110 meV. The dispersion of 
valence and conduction band is due to the hopping of 
electrons between neighboring molecules. We find that 
the Pc-derived eg conduction band and the empty big 
Cu-band are crossing but that there is very little hy- 
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FIG. 6: Band-structure, density of states and projected den- 
sity of states for (K2CuPc)2, calculated with SGGA-I-U. The 
band energies are plotted along lines in the first Brillouin zone, 
connecting the points in the Brillouin zone as before. 

bridization. 

Apparently the eg bands are also split, giving rise to 
a band gap. This is related to the fact that, as we al- 
ready noted above, in a free CuPc molecule the states 
eg symmetry are two-fold degenerate 0. In the K2CUPC 
crystal this degeneracy is lifted due to the ionic and cova- 
lent molecular crystal fields. The 110 meV band splitting 
between the Pc-derived valence and conduction band, 
however, is only a fraction of the factual charge gap of 
K2CUPC. As in the pristine system we have to take the 
Coulomb interactions between electrons on the Pc-ring 
{Upc and its screening) into account in order to find the 
correct gap value. After doing so we find A = 1.4 ± 0.1 
eV. 

In order to check how robust these results are, we have 
performed the same calculations and analysis for an iden- 
tical crystal structure, but with all the potassium ions 
removed. We refer to this as the KqCuPc structure. In 
this case we find a very similar band-structure to the 
one of K2CUPC, only with the Fermi-level shifted (Fig. 
01. Apart from this shift, the only appreciable change in 
band-structure is that doping induces a small shift of 150 
meV to higher energy of the almost dispersionless Cu big 
band. From this we conclude that potassium intercala- 
tion gives rise to a simple band filling of the CuPc-related 
electronic states. 

The rigid band behavior is important when we con- 
sider further doping of electrons in K2.75CUPC as in this 
compound we have for each CuPc molecule an extra of 
0.75 electrons, compared to K2CUPC. We calculated that 
the extra 0.75 electrons shift up the Fermi-level in Fig. 
Elby 0.154 eV. This implies that the doped carriers start 
filling the eg Pc-derived conduction band and the Cu- 
derived big band stays empty: the occupation of the big 
orbital turns out to be 0.0152 electron per CuPc. From 
this we conclude that electronic transport properties of 
(K2.75CUPC) are fully determined by the Pc-derived eg 
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FIG. 7: Band-structure, density of states and projected den- 
sity of states for (KoCuPc)2, calculated with SGGA-I-U. 



width of 0.32 eV. (K2CuPc) is a magnetic insulator as 
well, with a gap of 1.4 ±0.1 eV that is almost completely 
due to molecular Coulomb interactions. Our most im- 
portant result concerns the bandwidth of the conduc- 
tion band of (K2CuPc). It is 0.56 eV, which is large 
for a molecular solid. As, due to the rigid band behav- 
ior, in (K2.75CUPC) the additional 0.75 electron per CuPc 
molecule fill the eg Pc-derived conduction band, this wide 
and partially filled band entirely determines electronic 
transport. On the basis of these band structure calcu- 
lations one expects that (K2.75CUPC) is a metal. Note 
that in molecular crystal polaronic effects and disorder, 
not included in our analysis, tend to localize charge car- 
riers, especially in the low doping regime. However in 
(K2.75CUPC) the large amount of doped charge carriers 
enables an effective screening of disorder, and the unusu- 
ally large bandwidth makes polaron formation less likely, 
so that delocalization of the carriers is strongly promoted. 



states. 

Conclusions We find that pristine CuPc is a magnetic 
insulator with strong electron correlations, a charge gap 
of 2.5 ± 0.1 eV and a Pc-derived valence band with a 



Acknowledgments We thank Monica Craciun, Serena 
Margadonna and Alberto Morpurgo for stimulating dis- 
cussions and for sharing experimental data with us prior 
to publication. 



[1] T. Ishiguro, K. Yama ji, and G. Saito, Organic Super- 
conductors (Springer; 2nd edition, December 1998). 

[2] L. Forro and L. Mihaly, Rep. Prog. Phys. 64, 649 (2001). 

[3] M.F. Craciun, S. Rogge, D.A. Wismeijer, M.J.L. den 
Boer, T.M. Klapwijk, and A.F. Morpurgo, AIP Confer- 
ence Proceedings 696, 489 (2003); M.F. Craciun et al., 
cond-niat/0401036 , and to be published. 

[4] E. Tosatti, M. Fabrizio, J. Tobik, G. E. Santoro, Phys. 
Rev. Lett 93, 117002 (2004). 

[5] A. Rosa and E.J. Berends, Inorg. Chem. 33 584 (1994). 

[6] M.-S. Liao and S. Scheiner, J. Chem. Phys. 114 9780 
(2001). 

[7] G. Kresse and J. FurthmuUer, Phys. Rev. B 54, 11 169 
(1996); G. Kresse and J. FirthmuUer, Comput. Mater. 
Sci. 6, 15 (1996). 

[8] C. J. Brown, J. Chem. Soc. A 2488 (1968). 

[9] A. I. Lichtenstein, V. I. Anisonow, J. Zaanen, Phys. Rev. 



B 52 R5467 (1995) 
[10] M. W. Schmidt, K. K. Baldridge, J. A. Boatz, S.T. El- 
bert, M. S. Gordon, J. H. Jensen, S. Koseki, N. Mat- 
sunaga, K. A. Nguyen, S. J. Su, T. L. Windus, M. 
Dupuis, J. A. Montgomery, J. Comput. Chem 14, 1347 
(1983). 

[11] M.B.J. Meinders, J. van den Brink, J. Lorenzana and 

G.A. Sawatzky, Phys. Rev. B. 52, 2484 (1995). 
[12] J. van den Brink, M.B.J. Meinders, J. Lorenzana, R. Eder 

and G.A. Sawatzky, Phys. Rev. Lett. 75, 4658 (1995). 
[13] G. Brocks, J. van den Brink, A. Morpurgo, Phys. Rev. 

Lett. 93, 146405 (2004). 
[14] I.e. HiU, A. Kahn, Z.G. Soos, R.A. Pascal, Chem. Phys. 

Lett. 327, 181 (2000). 
[15] T. Schwieger, H. Peisert, M.S. Golden, M. Knupfer and 

J. Fink, Phys. Rev. B 66, 155207 (2002). 



